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Abstract 

The time-odd triaxial relativistic mean held approach is developed and applied to the investi- 
gation of the ground-state properties of light odd- mass nuclei near the double-closed shells. The 
nuclear magnetic moments including the isoscalar and isovector ones are calculated and good agree- 
ment with Schmidt values is obtained. Taking 17 F as an example, the splitting of the single particle 
levels ( around 0.7 MeV near the Fermi level ), the nuclear current, the core polarizations, and 
the nuclear magnetic potential, i.e., the spatial part of the vector potential, due to the violation of 
the time reversal invariance are investigated in detail. 
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I. INTRODUCTION 



Nuclear magnetic moment is one of the most important physics observables. It provides a 
highly sensitive probe of the single particle structure, serves as a stringent test of the nuclear 
models, and has attracted the attention of nuclear physicists since the early days jl|, 0,0]. Up 
to now, the static magnetic dipole moments of ground states and excited states of atomic 
nuclei throughout the periodic table have already been measured with several methods 
With the development of the radioactive ion beams (RIB) technique, it is now even possible 
to measure the nuclear magnetic moments of many short-lived nuclei near the proton and 
neutron drip lines with high precision 

Lots of efforts have been made to describe the nuclear magnetic moment non- 
relativistically or relativistically. In the non-relativistic theory, it has been pointed out 
early that the single-particle state of the shell model can couple to more complicated 2p-lh 
configurations and there are mesons exchange corrections caused by the nuclear medium 
effect 0. With modified g-factor or configuration mixing effect, non-relativistic approaches 



have turned out to be successful in reproducing data [2j. Therefore the single particle picture 
in the mean field approach may not be expected to describe the magnetic moment well . 
For LS closed shell nuclei plus or minus one nucleon, however, there is no spin-orbit partners 
on both sides of the Fermi surface and the magnetic moment operator can not couple to 
magnetic resonance jjj. Furthermore, the contributions from the pion-exchange current as 
well as others processes, e.g. A-hole excitation, to iso-scalar current have turned to be very 
small in the first order approximation QQ. Therefore, the iso-scalar magnetic moments 
of the light nuclei near double-closed shells should be described well in single-particle picture. 

In relativistic approach, the single particle wave function of nucleon is described by Dirac 
spinor with the large and small components, resulting in the relativistic effect in nuclear 
current and magnetic moment. Furthermore, as the spin-orbit coupling appears naturally 
in the relativistic approach, it is expected that the anomalous magnetic moment, which is 
related to the spin of nucleons, can be described well in relativistic approach with the free 
g-factor for the nucleon. 



During the 
success Il2. 13 



ast two decades relativistic mean field (RMF) theory has achieved great 



14, 



151 ] in describing many nuclear phenomena for stable nuclei, exotic nuclei 



HQ, „ we,, as snpernova and neutron stars Q. The RMF theory incorporates man y 
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important relativistic effects from the beginning, such as a new saturation mechanism by the 
relativistic quenching of the attractive scalar field, the existence of anti-particle solutions, 
the Lorentz covariance and special relativity which make the RMF theory more appealing 
for studies of high- density and high temperature nuclear matter, the origin of the pseudospin 



symmetry 



nucleon spectrum 




as a relativistic symmetry and spin symmetry in the anti- 

etc. Recently, good agreements with existing data for the ground-state 
properties of over 7000 nuclei has been achieved in the RMF+BCS mo del Q. 

However, it was found out that a straightforward application of the single-particle rela- 
tivistic model, where only sigma and the time-component of the vector mesons were con- 
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sidered, can not reproduce the experimental magnetic moments 
is attributed to the small renormalized nucleon mass (M* ~ 0.6M) which results in the 

n 

enhancement of the relativistic effect on the Dirac current 31]. As an improvement, the 
vertex corrections have to be introduced to define effective single-particle currents in nuclei, 
e.g., the "back-flow" effect in the framework of the relativistic extension of Landau's Fermi- 



liquid theory |32| or the RPA type summation of p-h and p-n bubbles in relativistic Hartree 
approximation |10L l^. 

In the widely used RMF approaches, there are only the time-even fields which are essential 
to physical observable. In odd-A or odd-odd nuclei, however, the Dirac current due to the 
unpaired valence nucleon will lead to the time-odd component of vector fields, i.e., the 
nuclear magnetic potential. In other words, the time-odd fields will have the same effect as 
the "back-flow" effect or the RPA type summation of p-h and p-n bubbles and give rise to 
the core polarization which will modify the nuclear current, single-particle spin and angular 
momentum, giving the appropriate magnetic moments. In fact, the time-odd fields are very 
important for the description of the magnetic moments 3^|, and rotating nuclei (313, 
etc. 

It is thus essential to consider the spatial-component of vector fields in the framework of 
RMF theory to investigate the single-particle properties in odd-mass nuclei. With the time- 
odd nuclear magnetic potential in axial deformed RMF, the iso-scalar magnetic moment have 
been reproduced well 0,13]. Therefore, it is interesting to include the time-odd potential 
into triaxial case [36 1 and investigate the distribution of non- vanishing spatial-component 
of the uj vector meson field, nuclear current, magnetic potential and magnetic moments in 
nuclei with odd numbers of protons and/or neutrons. Taking 17 F as an example, the density 
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distribution, single-particle energy and its splitting due to the violation of the time reversal 
invariance, nuclear potential, orbital momentum, etc., will be investigated. The importance 
of the time-odd component in the description of the nuclear current and magnetic moment 
will be illustrated. 

The paper is arranged as the following. In SecJHl the time-odd triaxial RMF approach is 
presented in detail. The splitting due to the violation of the time reversal invariance in the 
single-particle energy will be estimated by reducing the Dirac equation with the time-odd 
nuclear magnetic potential to Schrodinger equation. The numerical details for solving the 
triaxial RMF equations expanded in three dimensional harmonic oscillator basis are given in 
Sec JlllI The results and discussions are presented in Sec JIVl and a short summary is given 
in SeclVl 



II. TIME-ODD TRIAXIAL RELATIVISTIC MEAN FIELD APPROACH 



The starting point of the RMF theory is the standard effective Lagrangian density con- 
structed with the degrees of freedom associated with the nucleon field (•0), ty\ 
meson fields (a and u), the isovector meson field (p) and the photon field (A) |l2, 



13. 



C 



i^dp — M — g a a — gjtUp - g P l^r ■ p M - e'fA, 



1-T 3 



— - F F^ v 



(1) 



in which the field tensors for the vector mesons and the photon are respectively defined as, 





= d^u u 


- d u u^, 




= d^p u 




F 


= d»A v 


-d v A„ 



(2) 



The Lagrangian ((!} includes respectively the non-linear self-coupling of the cr-meson and 
the cu-meson. In this paper the arrows are used to indicate vectors in isospin space and the 
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bold types for the space vectors. The parameters include the effective meson masses, the 
meson-nucleon couplings and the non-linear self-coupling of the mesons. Because of charge 
conservation, only the 3rd-component of the iso- vector p-meson contributes (for brevity the 
superscript '3' is dropped out in the following). 

From the Lagrangian in Eq.([T|). the equations of motion for nucleon is 



{a . [_ zV - V(r)] + (3M* + Vb(r)}^(r) = e^(r) 
where M*(r) is defined as 

M*(r) =M + S(r), 
with the attractive scalar potential S(r), 



(3) 



(4) 



S(r) = g a a(r), (5) 
and the usual repulsive vector potential, i.e., the time component of the vector potential, 

1 — r 3 

^o(r) = &M)(r) + fiV^PoO) + e — - — A (r). (6) 
The time-odd nuclear magnetic potential, i.e., the spatial components of the vector potential, 
V(r), 

V(r) = g u u(r) + g p r 3 p(r) + e^pA(r). (7) 
is due to the spatial component of the vector fields. Compared with u>(r) field, p(r) and A(r) 

fields turned out to be small in V(r) so that they were often neglected for light nuclei 
The Klein-Gordon equations for the mesons and the electromagnetic fields are 



-V 2 + m?)C(r) = S c (r) 



(8) 



where, 



S ( (r 



9aps(r) - g2° 2 (r) - £3^ 3 0), ( = a 



9w3 bO) - c 3 uj u uj u ujV(r 



(9) 



9p3r(* 



e#(r), 



t = A» 



with, 

Ps(r) = Y, n i$i( r )' l Pi( r ), 

i 

f B ( v ) = E n ^i( r )7^i(r), 

* - (10) 
Jr( t ) = E^Oh^s^Mr), 

i 

l _ T .g 

Jc( r ) = E^( r )7 M — 5— V'tW. 

The sums are taken over the particle states only, i.e. the contributions from the negative- 
energy states are neglected (no-sea approximation), and rii is the occupancy probability of 
state ipi{r), which satisfies the normalization condition 

! rf 3 r^n^j(r)^(r) = Z (or N), (11) 

J i 

where Z (N) is the proton (neutron) number. If the pair correlation is neglected, the occu- 
pancy rii is one or zero for the states below or above the Fermi surface. 

The total energy of the system, including the spatial-component of vector fields, is, 



-Etotal = -Epart. + + E u + E p + E y + E cm + -Ep a i r , (12) 



where, 



E 



part. 



i=l 



E a 


= ~\) 


d 3 r 


g a p s (r)a(r) H 


~ ^2^ 3 (r) 


1 4 

+ 2 9sa 


(r) , 


E„ 


--W 


d 3 r 




- g u <*>(r) 


•jBW- 




E P 


=-\\ 


d 3 rg pPo (r)j° R (r), 








E~ 


e 2 
8tt 


1 dh 


A)(r)i e °(r), 









(13) 



with Ei = €i — M and E pai i V is due to the pairing correlation. 

The center-of-mass correction to the binding energy from projection-after-variation in 
first-order approximation is given by 
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Er 



2mA 



(P 2 ). 

\ cm/ 



(14) 



The root-mean-square (rms) radius is defined as, 



2\l/2 



d rr p(r) 
d 3 rp(r) 



1/2 



(15) 



The quadrupole moments Q^f and f° r neutron and proton are respectively calculated 
by 



5 

V20 



167T 



(2s 3 -j/ 2 -* 2 ), 



n,p 
22 



15 

32^ 



(^ 2 -y 2 ). 



(16) 



The deformation parameters f3 and 7 can be obtained respectively by the corresponding 
quadrupole moments for neutron, proton and matter (e.g., Q20 = Q20 + Q20) as, 



/5 = 77^\/Qio + 2gi 2 , T-tan- 1 ^-^ 



3A.R 2 







(17) 



20 



with Q c 



16tt 



Q20 + 2< 522 and Rq = 1.2 A 1 / 3 fm. 



Before any sophisticated numerical calculation, it is interesting to investigate the time- 
odd potential in the Dirac equation qualitatively by reducing the Dirac equation to more 
familiar Schrodinger-like form and estimate the splitting of the time-reversal orbits due to the 
time-odd potential g u (jJ. The Dirac equation Eq.fjSJ) can be reduced to the Schrodinger-like 
equation for the large component as, 



[tr ■ (p - g^) 



1 



e + M + S-V 



a ■ (p - g u u) + V + S}f = (e - M)f. 



(18) 



The time-odd part of the Hamiltonian can be obtained as the following: 



H. 



odd 



9u 



E + 2M + S - V 



[a- ■ u)(tr ■ p) - a- ■ p 



9u 



E + 2M + S — Vq 



a ■ u, (19) 
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Introducing uj = x r ; which defines a nuclear magnetic field B n = Vxw, the time-odd 
part of the Hamiltonian can be rewritten as: 

H^aa= — (L + 2S)-B n , (20) 

odd E + 2M + S-V K J ' K J 

where the assumption uj x p = is used and the higher order term . 



(E + 2M + S- Vo) 2 ' 
is neglected. 

It indicates that the spatial-component of vector meson field uj will lead to a time-odd 

term in the Schrodinger-like equation and result in the energy splitting. The Kramer's 

degeneracy in odd-mass nuclei is thus removed by the effective intrinsic nuclear magnetic 

dipole moment, — — — — (L + 2S), coupled to the nuclear magnetic field B n . The 
E + 2M -\- o — Vq 

magnetic field B n , as shown in the following, is around B n ~ 5x 10 -4 fm -2 in 17 F, which will 
lead to the energy splitting between two time reversal conjugate states, e.g., ~0.4 MeV for 
lsi/2- Of course, the more realistic splitting will depend on the detailed information on the 
nuclear magnetic field B n , the scalar and vector potentials and the density distribution of the 
particular orbits, which can be obtained by solving the Dirac equation © self-consistently. 



III. NUMERICAL DETAILS 

For the triaxial deformed nucleus, the Dirac spinor can be expanded by using the three- 
dimensional harmonic oscillator wave function $ Q (r, s) and $^(r, s) in Cartesian coordinates: 



$o,(r, s) = <f) nx (x)(f) n V (y)(f>n z ( z )—K I witn °L= \n x ,n y ,n z ,m s ). 

V2 \ / 

i h y ( 1 

$a(r,s) = (phA^hy^nXz)—, = (-l) n * +n y +1 I with a = \n x , n y , n z , m s ), 

V2 \ (-!)«■ 



(21) 



where, the phase factor i Uy is chosen in order to have a real matrix elements for the Dirac 
equation 38]. The time reversal conjugate states $o(r, s) and $^(r, s) satisfy, 

f $«(r, s) = -$ a (r, s), f$ = (r, s) = $„(r, s). (22) 

with T = —ia y K. 
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The normalized oscillator function <fi nk {k) in the k direction (k denotes x, y, or z) is, 



^W = %^(^)exp[-i(A) 2] 

V»fc °k ^ Ok 



(23) 



with the normalization factor N nk = (^/Tr2 nk n}.\) 1//2 . The Hermite polynomial H n (^) is 
given by 

Hrtt) = (-lTe e ^e-Z 2 . (24) 
The frequency in fc-direction Uk can be written in terms of the deformation parameters (3 
and 7 as, 

/ 5 2kn 
hw k = hhj exp[J— /5cos(7 — )], (25) 



with oscillator length b k = a/ h/ Mu>k and M the mass of the nucleon. The oscillator fre- 
quency is Jtwq = 41y4 -1 / 3 MeV and the corresponding oscillator length b = a/ H/Muq. 
The Dirac spinor for nucleon has the form 



| /(r ' S) 
ig{r,s) 



(26) 



where Xiif) ls the isospin part. In odd-A nuclei, as the time reversal invariance is broken 
by the unpaired valence nucleon, Eq.([26j) can be written as the linear combination of time 
reversal conjugate basis 



f( r > s ) = E4I«) + T,h\oc), 

a a 



(27) 



However, one can define the time reversal conjugate states ipj = Tipj and ih = —Tipj as, 



ipj(r,t,s) 



a 



J 



Xu (*), 



Xt 7 (t), 



(28) 



so that the Dirac equation for the nucleons can be solved separately in subspaces {ijij} or 
{^j}, i-e., 



and 




More details can be found in Appendix. 

Using the relations TcrT^ 1 = —er and TVT -1 = V, one can have 

(a\a ■ iV(r)|a) = -{a\a ■ iV(r)\a), (31) 

which will lead to different Hamiltonian matrix elements in subspaces {ipj} and {ipj}, as 

shown in Eqs. (J29J) and (JBTij) . Therefore the time-odd nuclear magnetic potential in Eq.(J7J) 
will result in the violation of the time reversal invariance and the splitting between the time 
reversal conjugate states Ej and Ej. 

The meson fields £(r) are also expanded in three dimensional harmonic oscillator basis 
with the same deformation parameters (3 and 7, 

,n y ,n z ( t ) n x {. X ) ( t ) n y {y) ( t ) nX Z )- (32) 

f^x i^Ly i^z 

In order to simplify the calculations and to avoid additional parameters, the oscillator length 

for mesons bs is smaller by a factor of \/2 than that bo for the nucleons. The Klein-Gordon 
equations (jHJ) for the meson fields become the inhomogeneous sets of linear equations, 

N B , 

M n'Wl<,n x n y n z C^n y ,n z = J dxdydZ(j) n n{x)^ n n{y)^ n >>{z)S^X,y,z), (33) 

with the matrix elements, 
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M. » a a = \\fruXru- — 1)5 » n + \firu- + l)(?v + 2)5 » , 15 » 5 » 

* n x n y n z ,n x n y n z 7)2 L V '"xv'x /"ri^n^— 2 1 V V ' v x 1 -vvx i "/"n^ na;+2J n y n y n z n z 

1 11 11 11 

(34) 

The source terms in Eq.fjlOjl can be obtained from the nucleon wave functions in two sub- 
spaces. For the Coulomb field, due to its long range character, the standard Green function 
method is usedj^. 

The expansion for the spinors and meson fields has to be truncated to a fixed number of 
basis states and the cut-off is chosen in such a way that the convergence is achieved for the 
binding energies as well as the deformation parameters (3 and 7. Here in the calculation, up to 
nf = 12 and rib = 10 are adopted, which gives an error less than 0.1% for the binding energy 
as demonstrated in Ref. js^J, and the effective interaction parameter sets PK1 |4lJ has been 



used. As shown already in Ref. [36| , the binding energy and the deformations of the nucleus 
are independent on the deformation of the basis in Eq. (25). Therefore for convenience, the 
deformation of the basis in Eq. (25) are chosen as (3 = 0, i.e. , = b . The center-of-mass 
correction is considered microscopically as in Eq. (jl4j) . The expanded oscillator shells for the 
small component g is carried out with n g = n/ + 1. In the calculation, the mirror symmetry 
with respect to the xy, xz, and yz planes for the densities and rotational symmetry for the 
currents are assumed to avoid the complex coefficients in the basis, which is appropriate as 
illustrated in Ref. (35), [3] . 



IV. RESULTS AND DISCUSSION 

Using Woods-Saxon potentials as initial potentials in Eqs.(jSJ) and © in the Dirac equa- 
tion©, and one tenth of the time-component of the vector meson field as the initial spatial- 
component, the above Dirac equation Q and Klein-Gordon equations (jHJ) are solved self- 
consistently by iteration. 



11 



Taking 17 F as an example, the importance for taking into account self-consistently the 
spatial-component of the uj field due to the violation of time reversal symmetry in odd-mass 
nuclei has been demonstrated by examining the nuclear magnetic potential, the nuclear 
magnetic field, the ground-state properties including the binding energy, rms radii and de- 
formation f3 and 7, the single-particle energy and the splitting of the time reversal conjugate 
states, the density and current distribution, etc. 

Finally, the nuclear magnetic moments for the light double-closed shells nuclei plus or 
minus one nucleon, including the Dirac, anomalous, iso-scalar, and iso-vector magnetic mo- 
ments, will be investigated and compared with data and former calculation results available. 

A. Nuclear magnetic potential 




FIG. 1: (color online) The distribution of the time-component of the vector potential and the 
nuclear magnetic potential, i.e., the spatial components of the vector potential, versus x axis for 
y = z = 1.29 fm in 17 F. 

There is no time-odd magnetic potential in the ground-state of even-even nucleus because 
of the time-reversal symmetry. For odd-odd or odd A nucleus, the unpaired valence nucleon 
will give non-vanishing contribution to the nuclear current which provides the time-odd 
magnetic potentials. The time-component of the vector potential and the nuclear magnetic 
potential ( the spatial components of the vector potential) in 17 F are given in Fig^ as 
functions of x for y=z=1.29 fm. Each component of the nuclear magnetic potential has 
a peak value around 1.0 ~ 1.5 MeV. Compared with the time-component of the vector 
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potential Vo, the nuclear magnetic potential is two order smaller in magnitude. However, 
as the nucleon moves in a potential ~ —50 MeV given by the cancelation of the attractive 
scalar potential S(~ —400 MeV) and the repulsive time- component of the vector potential 
Vo(— +350 MeV), the magnetic potential will play important role in single-particle properties 
as will be shown later. Furthermore, Fig^ shows that the magnetic potential is comparable 
in magnitude with the iso-vector vector potential provided by the time-component of p 
meson field. 

With the presence of the spatial-component of vector field, the magnetic field B n can 
be defined as B n = V x uj. The distribution of nuclear magnetic field B n versus x for 
y=z=1.29 fm in 17 F are shown in FigJ21 It is found that B n changes rapidly with x and its 
x-component is larger than the other two components. The value used in the estimation of 
the energy splitting between time reversal conjugate states is around one third of the peak 
value for the x-component of B n for y=z=1.29 fm in 17 F. 




FIG. 2: (color online) The nuclear magnetic field B n (= V x u) versus x axis for y=z=1.29 fm in 
17 F. 

Taking into account the time-odd potentials self-consistently will lead to differences 
around 3 MeV for E a and E u in Eq. (jl3|) . Moreover, the nuclear magnetic potential will 
lower the potential (Vo + S)(r) in 17 F as shown in Fig|Hl However, these differences will 
cancel with each other. Therefore the time-odd magnetic potentials have small influence 
on the binding energy, rms radii and quadrupole moment Qq for 17 F. The binding energy 
per nucleon, rms radii and quadrupole moment for 17 F in triaxial RMF theory with (with- 
out) the time-odd magnetic potentials are respectively 7.600(7.589) MeV, 2.63(2.63) fm and 
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9.78(9.86) fm 2 . The deformation parameters j3 and 7 for 17 F are respectively 0.08 and 57.14° 
after taking into account the time-odd magnetic potentials self-consistently. The (3 value is 
similar as the previous axial deformed RMF calculation, i.e., — 0.086 (42 1 . 
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FIG. 3: (color online) The potential (Vo + S)(r) for proton and neutron in 17 F with (solid line) 
or without (dot-dashed line) the time-odd magnetic potentials self-consistently in triaxial RMF 
theory. 

B. Splitting of the time-reversal conjugate single-particle states 

The violation of the time reversal invariance will result in the core polarization and lead 
to the splitting of the the time-reversal conjugate states. Moreover, Because of deformation, 
the single-particle angular momenta are no longer good quantum numbers. The energies 
and expectation values of the angular momentum for single-particle time reversal conjugate 
states ipj and ipj in 17 F are respectively given in Table |1J The expectation values of the 
orbital angular momentum for the large ((L)) and small component ((L')) of the single- 
particle state are close to their spherical values due to the small (3. As 7 is close to 60°, 
L x + S x is an approximate good quantum number and close to its axial deformed value, from 
which we can determine the approximate quantum number Ijm for the single-particle states 
as given m Table [Q in this case. One also has (L f ) ~ (L) ± 1 for the state J — L ± ^- As 
expected from the estimation in Secjnl the energy splitting for single-particle time reversal 
conjugate states in 17 F range from 0.04 to 0.7 MeV and larger splitting occurs for states 
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with larger expectation value (L x + S x ). 

TABLE I: The energies and expectation values of the angular momentum for single-particle time 
reversal conjugate states ipj and if?j in 17 F, where L and V are the orbital angular momenta of 
large and small components respectively, L x and S x are the x-component of the orbital angular 



momenta and spin, and Ae = £j — Ej. 



proton 


lj\m\ 










e (MeV) L (ft) L' (ft) L x (ft) 5 X (ft) 


e (MeV) L (ft) V (ft) L x (ft) S x (ft) 


Ae (MeV) 


S l l 
2 2 


-37.55 


0.01 


1.01 


0.00 


-0.50 


-37.92 


0.01 


1.01 


0.00 


0.50 


0.37 


P3 3 
2 2 


-17.93 


1.01 


1.93 


-1.00 


-0.50 


-18.60 


1.01 


1.92 


1.00 


0.50 


0.67 


P3 1 

2 2 


-17.92 


1.01 


1.97 


-0.41 


-0.10 


-17.88 


1.01 


1.97 


0.40 


0.10 


-0.04 


Pi 1 

2 2 


-11.60 


1.03 


0.28 


-0.59 


0.10 


-11.83 


1.03 


0.28 


0.60 


-0.10 


0.23 


ds 5 

2 2 


-1.14 


1.98 


2.85 


-1.98 


-0.50 


-1.83 


1.92 


2.83 


2.00 


0.50 


0.69 


C?5 3 
2 2 


-1.10 


2.02 


2.91 


-0.95 


-0.20 


-1.27 


2.03 


2.94 


1.22 


0.25 


0.17 


G?5 1 

■>. 9. 


-1.01 


2.01 


2.93 


-0.39 


-0.10 


-1.08 


2.02 


2.98 


0.12 


0.03 


0.07 


neutron 












e (MeV) L (ft) L' (ft) L x (ft) S x (ft) 


e (MeV) L (ft) V (ft) L x (ft) S x (ft) 


Ae (MeV) 


S 1 1 
2 2 


-43.72 


0.01 


1.01 


0.00 


-0.50 


-44.10 


0.01 


1.01 


0.00 


0.50 


0.38 


P3 3 
2 2 


-24.26 


1.01 


1.93 


-1.00 


-0.50 


-24.96 


1.01 


1.93 


1.00 


0.50 


0.70 


P3 1 
2 2 


-23.73 


1.01 


1.95 


-0.48 


-0.02 


-23.69 


1.01 


1.96 


0.47 


0.03 


-0.04 


pi 1 

2 2 


-17.40 


1.02 


0.33 


-0.52 


0.02 


-17.65 


1.02 


0.33 


0.53 


-0.03 


0.25 


C?5 5 
2 2 


-6.93 


1.93 


2.81 


-2.00 


-0.50 


-7.69 


1.93 


2.82 


2.00 


0.50 


0.76 


C?5 3 
2 2 


-6.35 


2.00 


2.91 


-1.23 


-0.21 


-6.54 


2.00 


2.93 


1.27 


0.22 


0.19 


G?5 1 

2 2 


-6.15 


1.97 


2.90 


-0.38 


-0.11 


-6.22 


1.98 


2.91 


0.34 


0.10 


0.07 



In order to investigate the difference between the energy splitting due to the violation of 
time reversal invariance from valence neutron and proton, the proton and neutron single- 
particle levels for time reversal conjugate states in 17 and 17 F are shown in Fig0] As a 
reference, those in 16 as well as in 17 and 17 F with time-even potentials only are also 
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given since 17 or 17 F is nothing but one neutron or proton on top of the 16 core. 

There are the following features: 1) The valence proton (neutron) in 17 F ( 17 0) will lower 
the neutron (proton) single particle levels more noticeably than the corresponding proton 
(neutron) single particle levels; 2) The energy splitting due to the violation of time reversal 
invariance from valence proton (neutron) in 17 F ( 17 0) is distinguishable in Fig0] which 
denotes the necessity to take into account the time-odd magnetic potentials. In Ref. 42], 
the influence of the time-odd magnetic potentials on single-particle energies and the related 
quantities such as single-particle separate energy has been investigated. 
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FIG. 4: (color online) The proton and neutron single-particle levels for time reversal conjugate 
states in 17 and F. Those in 16 as well as in 17 and 1 F without time-odd potentials are also 
given as a reference. The levels with negative angular momentum projection are given in blue (the 
3rd and 6th columns) and the levels with positive angular momentum projection are given in red 
(the 4th and 7th columns) . 

The density distribution of the unpaired valence nucleon for 17 and 17 F as well as the 
proton, neutron and matter for 16 0, 17 0, and 17 F in yz, zx, and xy plane are plotted in 
FigsEl and El respectively, in which the other axis x, y, and z is integrated respectively. It 
is interesting to observe that although there are noticeable differences in the energy of the 
unpaired valence nucleon, the density distributions of the valence proton in 17 F and neutron 
in 17 are quite similar to each other and they are corresponding to lc/5/2 orbits. In Fig® 
the density distributions for proton, neutron and matter in 17 and 17 F as well as 16 are 
shown. 
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FIG. 5: (color online)The density distribution of the unpaired valence nucleon in yz, zx, and xy 
plane for 17 F and while the other axis x, y, and z is integrated respectively. 
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FIG. 6: (color online) The proton, neutron and matter density distribution in yz, zx, and xy plane 
for 16 O, 17 O, and F, while the other axis x, y, and z is integrated respectively. 

C. Current and nuclear magnetic moments 

Nuclear magnetic moment is related to the effective electromagnetic current operator [3, 

Q 
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FIG. 7: (color online) Upper: The Dirac current for proton in yz plane at x=0.26 fm in 17 F. The 
direction and length of the arrows respectively represent the orientation and magnitude of current. 
Lower: The azimuthal component of Dirac current for proton as a function of z with y=2.32 fm 
and x=0.26, 0.77, 1.29, 1.80, 2.32 fm respectively in 17 F. 

P{x) = Q$(x)^{x) + ^d^(x)a^ij(x)}, (35) 

where the field operators are in the Heisenberg representation with Q = |(1 — r 3 ), M the 

nucleon mass, a^ u = f [7^,7^], and k the free anomalous gyromagnetic ratio of the nucleon: 
kP = 1.793 and K n = —1.913. The spatial-component of the current operator is given by, 

j(r) = Q^ (r)mKr) + _^_ v x ^t (r)/?s ^ (r)]; ( 36 ) 

where the first term gives the Dirac current joM — ip^{r)aip(r), which can be decomposed 

into a orbital (convection) current and a spin current according to Gordon identity, and the 
second term in Eq. ()36|) is the so-called anomalous current. 

The proton Dirac current in yz plane at x=0.26 fm in 17 F is given in the upper panel in 
Fig|7|in which the direction and length of the arrows respectively represent the orientation 
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and magnitude of current. It is clearly seen that the Dirac current peaks at the nuclear 
surface. In order to investigate how the Dirac current change with x, the azimuthal com- 
ponent of Dirac current for proton is given for x = 0.26, 0.77, 1.29, 1.80, 2.32 fm respectively 
as a function of z with y = 2.32 fm in 17 F. It is found that the Dirac current decreases with 
x. Compared with the density distribution for the valence nucleon in Fig. |3J it is obvious 
that the current is mainly from the valence nucleon and moves around the nuclear surface 
in yz plane. In order to see the core polarization effect, the nuclear magnetic moment in the 
following is more suitable than the current. 

In the mean field approximation, the nuclear magnetic moment can be calculated from 
the expectation value of the current operator in Eq. (|3l)|) in ground state, 



/x = 2M J rf 3 r^[r x {g.s.\j(r)\g.s.)] 

A 



Mc 2 

-— QVl(r) r x a M r ) + M^~(r)/3£^(r) 



(37) 



Here, h and c are added in order to give the magnetic moment in units of nuclear magneton. 
The magnetic moment can be divided as the Dirac magnetic moment, 

A Mr 2 C 

= E 1" / d 3 rQ^(r)r x c^(r), (38) 

2 = 1 

and the anomalous magnetic moment, 

Ma = K ! d 3 rfi A (r) = « / rf 3 r^(r)/3S^(r). (39) 

i=X J 

The anomalous magnetic moment (proton) density [ia{?) distribution in x=0.26 fm plane 

and the distribution of /U^(r) (proton) along z in y=0.26 fm plane, with x=0.26, 0.77, 1.29, 
1.80, 2.32 fm respectively in 17 F are plotted in Fig|HJ It is found that the anomalous magnetic 
moment density distribution has the similar character to the Dirac current since they are 
both mainly dependent on the density distribution of the unpaired valence proton given by 
FigEl 

The magnetic moments in Eq. (j37j) as well as the Dirac and anomalous magnetic moments 
in Eqs. (J3~8"j) and (|3^j) for light nuclei near double-closed shells with A=15, 17, 39, 41 are 
respectively given in Tables HP and II I II 
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FIG. 8: (color online) Upper: The anomalous magnetic moment density [ia(t) f° r proton in yz 
plane at x=0.26 fm in 17 F. Lower: The anomalous magnetic moment density ha{ t ) f° r proton as 
a function of z with y=0.26 fm and x=0.26, 0.77, 1.29, 1.80, 2.32 fm respectively in 17 F. 



TABLE II: The magnetic moments of light nuclei near double-closed shells in units of nuclear 
magneton /ijy. The spherical and axial deformed calculations with NL1 are taken from Ref.jj]]. 



I 1 


15Q 


17 


39 Ca 


41 Ca 


15 N 


17 F 


39 K 


41 Sc 


Exp. 


0.72 


-1.89 


1.02 


-1.60 


-0.28 


4.72 


0.39 


5.43 


Tri.RMF 


0.57 


-2.00 


0.98 


-2.13 


-0.19 


4.893 


0.37 


6.04 


Axi.RMF 


0.65 


-2.03 


0.96 


-2.13 


-0.29 


4.99 


0.33 


6.07 


Sph.RMF 


0.66 


-1.91 


1.17 


-1.91 


-0.03 


5.05 


0.72 


6.32 


Schmidt 


0.64 


-1.91 


1.15 


-1.91 


-0.26 


4.79 


0.12 


5.79 



As shown in Table |Hj the magnetic moments given by time-odd axial deformed RMF 
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TABLE III: The Dirac and anomalous magnetic moments ( lid and \ia ) of light nuclei near closed 
shells in units of li-^- The spherical and axial deformed calculations with NL1 are taken from 
Ref.bl. 



fJ-D 


15/~> 

10Q 




jy Ca 


41 n 


1 5 AT 


1 7t 




41 o 

4i Sc 


Tri.RMF 


-0.11 


-0.13 


-0.16 


-0.28 


0.46 


3.15 


-0.64 


4.31 


Axi.RMF 


-0.13 


-0.13 


-0.30 


-0.22 


0.44 


3.21 


1.50 


4.29 


Sph.RMF 










0.59 


3.26 


1.81 


4.54 


Schmidt 










0.33 


3.00 


1.20 


4.00 


fJ-A 


15Q 


17 


39 Ca 


41 Ca 


15 N 


17 F 


39 K 


41 Sc 


Tri.RMF 


0.68 


-1.87 


1.13 


-1.85 


-0.64 


1.75 


1.00 


1.74 


Axi.RMF 


0.78 


-1.90 


1.26 


-1.87 


-0.73 


1.78 


-1.17 


1.79 


Sph.RMF 


0.66 


-1.91 


1.17 


-1.91 


-0.62 


1.79 


-1.09 


1.79 


Schmidt 


0.64 


-1.91 


1.15 


-1.91 


-0.60 


1.79 


-1.08 


1.79 



TABLE IV: Iso-scalar magnetic moments fJ,(A) = \p(Z, N) + fj,(Z + 1, JV — l)]/2 of light nuclei near 
closed shells in units of fx^. 



A Orbit EXP. Tri.RMF Axi.RMF [9] Sph.RMF[9] Schmidt 



15 lpx/2 


0.22 


0.19 


0.18 


0.32 


0.19 


17 ld 5/2 


1.41 


1.45 


1.48 


1.57 


1.44 


39 ldg/a 


0.71 


0.67 


0.64 


0.94 


0.64 


41 I/7/2 


1.92 


1.96 


1.97 


2.21 


1.94 



and triaxial deformed RMF calculations reproduce the data well in most cases and are close 
to each other due to the small 7 deformation for the nuclei investigated here. The relative 
large difference for 15 N in Table HP between the axial and triaxial RMF calculations is mainly 
from the anomalous part and it is due to the inclusion of 7 degree of freedom. 

The anomalous magnetic moments given by the triaxial RMF calculations are close to 
those of the spherical ones. While the difference in the total magnetic moment between 
the spherical and triaxial calculations is mainly from the Dirac magnetic moments due to 
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TABLE V: Iso- vector magnetic moments n(A) = \p(Z, N) — \i{Z + 1, N — l)]/2 of light nuclei near 
closed shells in units of /xn- 



A 


Exp. 


Tri.RMF Axi.RMF[9j 


Sph.RMF[9j Schmidt 


15 


0.501 


0.376 


0.470 


0.345 


0.451 


17 


-3.303 


-3.446 


-3.510 


-3.480 


-3.353 


39 


0.312 


0.305 


0.315 


0.225 


0.512 


41 


-3.513 


-4.086 


-4.10 


-4.115 


-3.853 



the core polarization. Without core polarization, the Dirac magnetic moment of neutron 
particle (hole) should vanish. However, as seen in Table the Dirac magnetic moment 
for 15 0, 17 0, 39 Ca, and 41 Ca contribute around 10% to the total nuclear magnetic moment 
which is the effect of the core polarization. The core polarization can be taken into account 
self-consistently in the axial and triaxial calculations but not in the spherical one. 

Although the Dirac and anomalous magnetic moments for 39 K given by the triaxial RMF 
calculation have different signs with the other calculations, the triaxial RMF calculation bet- 
ter reproduce the experimental magnetic moment. The triaxial deformed RMF calculation 
result of 39 K gives the deformation parameters f3 = 0.038 and 7 = 58.66°. 

The iso-scalar and iso-vector magnetic moments from the spherical, axial and triaxial 
deformed RMF calculations for light nuclei near the closed shells are given in Tables ITVl and 
IV1 in comparison with the data and Schmidt values. As mentioned before, the iso-scalar 
magnetic moments in deformed RMF calculations are close to the Schmidt values if the 
core polarization effect has been taken into account self-consistently. The spherical RMF 
calculations, however, have discrepancies around 10 ~ 30% with the data. For the iso-vector 
magnetic moments, the axial and triaxial deformed RMF calculations agree better with the 
data than the spherical RMF calculations and the Schmidt values, which demonstrates the 
importance of the core polarization. In fact, the agreement with the experimental iso-scalar 
and iso-vector magnetic moments can be improved also in the spherical RMF calculations 
after considering the core polarization by either the Landau-Migdal approach [3^] or the 
configuration mixing j^ . 
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V. SUMMARY 



The time-odd triaxial relativistic mean field approach has been developed and applied 
to the investigation of the ground-state properties of light odd-mass nuclei near double- 
closed shells. The splitting due to the violation of the time reversal invariance in the single- 
particle energy has been estimated by reducing the Dirac equation with the time-odd nuclear 
magnetic potential to Schrodinger equation. 

The Dirac equation with time-odd nuclear magnetic potential for nucleon and Klein- 
Gordon equations for meson fields have been solved self-consistently by expansion on the 
three dimensional harmonic oscillator basis with PK1. Taking 17 F as an example, the ground- 
state properties including the binding energy, deformation (3 and 7, the single-particle energy 
and the splitting of the time reversal conjugate states, density and current distribution, and 
magnetic moments, etc., have been examined. It is found that although the non- vanishing 
spatial-component of the uj field due to the violation of time reversal symmetry in odd-mass 
nuclei has small influence on the binding energy, root-mean-square radii and quadrupole 
moment, it will create a magnetic potential, change the nuclear wave function and result 
in the core polarization effect which plays important role on the single-particle properties 
and the magnetic moment. The nuclear magnetic moments for the light double-closed shells 
nuclei plus or minus one nucleon, including the Dirac, anomalous, iso-scalar, and iso-vector 
magnetic moments, have been investigated and good agreements with the data have been 
achieved. 

With the time-odd triaxial relativistic mean field approach developed here, it can be 
expected to describe well not only the light double-closed shells nuclei plus or minus one 
nucleon, but also the nuclei in the 7-soft mass region. Investigation along this line is in 
progress. 
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APPENDIX A: THE BASIS FOR SOLVING THE DIRAC EQUATION 

Using the eqs. ()26|) and (}2Tj) . the Dirac equation Q with time-odd magnetic potentials 
can be written in matrix form, 
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A a i a A a ' a $a'a B a i a \ ( f a ^ 
Aa'a Aa'a Ba'a $a'a fa 



ffA 

fa' 
9(3' 

\9-p> ) 



(Al) 



If the mirror symmetries with respect to three planes (xy), (yz) and (xz) are assumed for 
all scalar potentials and vector potentials j^sj , it can be shown that the matrix elements in 
script letters will vanish. For example, the invariance of the matrix element A a 'ai 



A a 'a 



1 _|_ (—l)K+* i x+ 1 



\n x +h y +l 



l y- ,l y 



-1 — 



d*r(j> n i{x)<l> n >{y)<l> nlx {z){y Q + M*)(f>n x (x)4>f. (y)<fih z (z) 



(A2) 



under the operation x — > —x will leads to A a ' a = 0. After similar operation for y and z, the 
matrix in Eq. ljAlj) can be written as, 



/ 


A , 


B a i a 








\ 


(fa\ 




(fa'\ 




p a 












9a 


= e 












Aa'a 


Ba'a 




h 




fa' 


V 








— Bp' a 


Cf3' a 


J 


\9a) 




\9p> J 



(A3) 



which indicates that the space {ipi} can be reduced into two subspace {ipj} and {^j in 
Eq.fZBD- 
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